SHEP 95-21 



Momentum Scale Expansion 

of 

Sharp Cutoff Flow Equations 



Tim R. Morris 

Physics Department 
University of Southampton 
Southampton, SOU IB J, UK 



Abstract 

We show how the exact renormalization group for the effective action with a sharp mo- 
mentum cutoff, may be organised by expanding one-particle irreducible parts in terms of 
homogeneous functions of momenta of integer degree (Taylor expansions not being pos- 
sible). A systematic series of approximations - the 0(p M ) approximations - result from 
discarding from these parts, all terms of higher than the M th degree. These approxi- 
mations preserve a field reparametrization invariance, ensuring that the field's anomalous 
dimension is unambiguously determined. The lowest order approximation coincides with 
the local potential approximation to the Wegner-Houghton equations. We discuss the 
practical difficulties with extending the approximation beyond 0(p ). 



hep-th/9508017 
August, 1995. 



1. Introduction. 

In ref.[|l], Wegner and Houghton formulated an exact renormalization group equation 
for the Wilsonian effective actionB Sa[<£>] where the associated momentum cutoff A is taken 
to be sharp: only momentum modes p satisfying p < A are kept in the functional integral. 
The flow equation corresponds to computing the induced change in S\ as A is lowered 
to A — SA, by integrating out the modes with A — 5A < p < A. In the same paper, 
the authors solved for the Wilson fixed point (i.e. massless case) of the large N limit 
of three dimensional 0(N) scalar field theory^, using the exact renormalization group 
language. This solution involved approximating the effective action by just a general 
effective potential V\(<p): 



S A [<p] = Jd D x ±(d^) 2 + V A &) ■ (1.1) 



In the large N limit this approximation effectively becomes exact. In ref.0, Nicoll, Chang 
and Stanley proposed to use such a local potential approximation even when there is 
no justification in terms of the smallness of 1/N. Such a local potential turns out to 
be one-particle irreducible and may equivalently be regarded as an approximation to the 
sharp cutoff flow equations for the Legendre effective action |5|,|l3|].i This approximation 



has proved to be very robust and has been employedS many times since [|]-|12|. Obviously, 
it is very important to be able to establish such approximations, because there are many 
situations in quantum field theory where no genuinely small parameter exists to control 



the approximation [|T9 1 . 



One of the main reasons for the present paper is to provide further theoretical justifi- 
cation for the local potential approximation by demonstrating that it may effectively be 
regarded simply as the start of a certain expansion of S\ in powers of momenta. This is by 
no means as trivial as it might appear at first sight. Quite apart from the fact that a direct 
expansion in small momenta of the one-particle reducible parts of S\ vanishes to all orders, 
a Taylor expansion in momentum components p^ cannot any way be implemented: the 
sharp cutoff in the flow equations induces non-analyticity at the origin of momentum space, 
reflecting the non-local behaviour in position space. Instead, the expansion is performed 



1 a.k.a. Hamiltonian in statistical mechanics language 

2 i.e. the large N limit iV-vector model, equivalent |2| to the exactly solved spherical model] 

3 a.k.a. course-grained Helmholtz free energy. This argument is reviewed later. 

4 and sometimes rediscovered. For derivations see e.g. refs.QMf^ll. 
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via the one-particle irreducible parts, and the best that one can achieve is an expansion in 
momentum scaie ||13|| p ~ y/p^Pfi , the M th approximation resulting from dropping all terms 
beyond 0{p M ). This raises the spectre that the lowest order approximation - 0(p ), in 
which the momentum scale dependence is completely discarded, nevertheless still involves 
Green functions with non-trivial dependence on the angles between their momenta. It 
turns out that this dependence cancels however, and it is in this way that we recover the 
local potential approximation PHT2|| . 

Our second main motivation is unfortunately more negative. We wish to provide the 



evidence for the claim in ref.[n]: that the momentum scale expansion of the sharp cut- 
off flow equations is not a practical technique, i.e. beyond O(p ), because only certain 
truncations (of the field dependence) are calculable in closed form, and this is not good 
enough|8|]. It is not good enough because, in contrast to a momentum scale expansion - 
which might be expected to workU since it corresponds to an expansion in (an appropriately 
defined) 'localness' of the effective Lagrangian, truncations of the field dependence have 
limited accuracy and reliability [Q. We want to emphasise here that our claim that there 
are practical difficulties in extending the momentum scale expansion beyond O(p ), merely 
amounts to an admission of our own lack of ability in extending it in a substantive way. It 
would certainly be worthwhile to improve on these attempts (presented in sect. 6), because 
a simple two-loop example - and general arguments (c.f. ref.|l3[] and appendix B), suggest 
that the sharp cutoff momentum scale expansion should yield the fastest convergence (i.e. 
faster than the smooth cutoff momentum expansions fl~4l , |l5| ) . 

Despite the increasing interest in utilising the exact renormalization group to develop 
non-perturbative approximations , there have so far been relatively few works concerned 
with going beyond the local potential approximation |TT|-[T8[1 . Part of the point of sect. 6 
is also pedagogical: we provide simple concrete illustrations of the discussions in earlier 
sections, but also reemphasise [ p!T|p^| - |T6| , |20[] in a simple context the important issue of field 
reparametrization invariance. It has been understood in the condensed matter literature for 
quite some time|2(],|16] that it is important^ to preserve a field reparametrization invariance, 
e.g. simply 

<p^<p/\ , (1.2) 



5 and certainly seems to so far [11] 

6 for linear renormalization group transformations such as are exclusively discussed in the 



present context |19|| 
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in the flow equations. For example, if field reparametrization invariance is broken by 
the approximation scheme, the field's anomalous dimension can no longer be determined 
uniquely - the value depends on at least one unphysical parameter. Approximations by 
expansion in momenta (equivalently derivative expansion in the smooth cutoff case) gener- 



ally break reparametrization invariance [|16| . To date, only two forms of cutoff function are 
known that allow a field reparametrization invariance after approximation by momentum 
expansion: the sharp cutoff studied here, and power-law smooth cutoffs |fH|. 

Our formulation at first closely follows that of ref.[0|, but we will organise things differ- 
ently. Therefore we feel it is preferable to give as much as possible a self-contained account, 



indicating the similarities or differences to ref. |T3j as they arise. In sect. 2, we rederive the 
sharp cutoff flow equations, taking due care of the limit implied, but also incorporate ex- 
plicitly the possibility of an anomalous dimension, and display the field reparametrization 
invariance. In sect. 3, the momentum scale expansion, and the approximations which it 
suggests, are defined - and we stress again the fact that Taylor expansions in momenta 
cannot be implemented. This means that the important property of locality needs to be 
recovered, and we do this in sect. 4; we also explain why, and how, the uniqueness of the 
expansion then follows. In sect. 5 we work out the approximation to 0(p ) and demonstrate 
that this coincides with the local potential approximation^. Sect. 6 addresses by example 
the new issues that arise in going beyond 0(p ), computing several examples, in particular 
we indicate why we seem from the practical point of view to be limited to truncations of 
the field dependence, and discuss - with reference to a simple truncation - the issue of 
field reparametrization invariance. Finally, in sect. 7 we present our summary, making 
some general observations on the momentum expansion presented here. 

Before embarking on the paper proper, we review below the arguments that lead to 
expressing the flow in terms of a Legendre effective action Ta^] with I.R. (infra-red) cutoff 
A. The reason that this, or an equivalent formulation, must be used is that the Polchinski 
(or equivalently 1 114|| the Wilson) effective action 5a. [<£>], for the effective theory with U.V. 
(ultraviolet) cutoff A, has a tree structure, composed of full propagators with I.R. cutoff 
A, and one-particle irreducible parts which are those generated by Ta [f] ■ This structure is 
already manifest in the graphical form of the Polchinski (equivalently [|nj Wilson|2!J, or 
in the limit of sharp cutoff, Wegner'sJIJ) equation for the vertices of S\, which is displayed 
in fig. 1. (These trees take account of the purely classical fluctuations that are integrated 
out as A is reduced, and would be there even in an effective action for the classical field 
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theory, while the one-particle irreducible parts summarise the integrated out quantum field 
fluctuations) . 




Fig.l. The Polchinski equation for the vertices of Sa- The vertices are drawn as open 
circles. In the limit of sharp cutoff, the black dot represents a delta-function restriction of 
the momentum q, in the propagator, to q = A. 



To resolve the limiting procedure implicit in taking a sharp cutoff, it is necessary to 
take this structure into account, which has the effect of reducing the equations to those 
(or closely equivalent to those) for T\[(p] (c.f. JTB|,pB|). Yet again, we must preserve 



the tree structure on taking a momentum expansion. The momentum expansion naively 
corresponds to Taylor expanding Sa in the scale of the external momenta (to the vertices 
of Sa), regarding this as small compared to the cutoff A. In the sharp cutoff limit, this 
would cause all tree terms with internal propagators^ to vanish however, since the internal 
propagator is furnished with a sharp I.R. cutoff A and the momentum flowing through 
this propagator is of the same scale as the external momenta by momentum conservation. 
Clearly this is too great a mutilation of the theory, in particular not only do all tree level 
corrections get discarded in the process, but all loop diagrams with more than one vertex 
get discarded also (since these arise from substituting the tree parts of Sa into the second 
term in fig. 1). Instead we apply a momentum expansion only to the one-particle irreducible 
vertices, equivalently to Ta- In these cases all internal propagators are integrated over 
momenta greater than A (as follows from integrating fig. 1 between the overall U.V. cutoff 
Aq and A), and the momentum expansion corresponds to expanding the external momenta, 
regarded as small compared to these internal momenta. (Actually this oversimplifies a 
little: the structure is hierarchical, with 'inner' loops containing momenta larger than those 



7 such as that of the first term in fig. 1 
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of the 'outer' loops and inner propagators expanded in the outer loops momenta |fL3! ). In 
this way no propagators are lost, the diagrammatic structure is respected, and we might 
reasonably hope that the momentum expansion leads to a convergent numerical! series. 
The model two-loop calculation presented in ref. |13J (and reviewed in appendix B) is 
encouraging in this respect since it results in a rapidly convergent numerical series. 

Finally we briefly mention an immediate consequence of the tree structure: Setting 
the external momenta for the vertices of S\ to zero, kills all the tree terms. Therefore 
the Wilson effective potential (supplied with U.V. cutoff A) coincides with the effective 
potential in the Legendre effective action, if the latter is computed with I.R. cutoff A [13]. 
This explains our statement in the first paragraph. 



2. The flow equations. 

Thus we take as our starting point the partition function defined as 

expVF A [J] = Jvip ex${-\V'C- 1 .ip-S A f i [ip\ + J.^} , (2.1) 

regularised by an overall U.V. momentum cutoff Aq. The notation is essentially the same 
as previously Hl3|Jl4j] , so two-point functions are often regarded as matrices in position or 
momentum (q) space, one-point functions as vectors, and contractions indicated by a dot. 
Momentum conserving 5-functions are factored out when appropriate. We work in D 
Euclidean dimensions with a single real scalar field (p. The definition differs from ref. ||13|| 
only in that here we include the kinetic term ^(9 M (^) 2 in S A f (and will do likewise in T\) 
so that Sif is the full bare action for the theory, while C~ 1 (q, A) = (l/8 £ (q, A) — 1) q 2 
is now an 'additive' I.R. cutoff. Here, 9 s (q, A) is a smooth regularisation of the Heaviside 
9 function, of width 2e, satisfying < 9 e (q, A) < 1 for all (positive) A and q, but with 
9 e (q, A) — > 9(q — A) as e — > 0. The additive form for the cutoff will make manifest the field 
reparametrization invariance: that is invariance of the flow equations under rescaling of 



the field (|1.2j). (This fact was already used in ref. [|T2|. The same effect would be obtained 
by redefining Fa i— > Ta — h f d D x (d^ip) 2 at the end of the calculation.) The reason for 
this can be understood as follows: For any finite e the invariance is broken by the induced 
change in the cutoff term C _1 i— > A 2 C _1 , but in the limit of sharp cutoff (i.e. e — > 0), 



8 It is not in fact an expansion in a small parameter because the external momenta p are 
integrated out over the range p < A. 
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A 2 C 1 has the same effect for any A, since it is either zero (for momenta q > A) or forces 
the integrand of (2.1) to vanish. From ( |2.1|) we have 



d_ 
dA 



W A [J] 



1 r SWa dC' 1 5Wa 

2 \1T' dA 



+ tr 



dC~ l 5 2 W A \ 
dA ' 5JSJ ) 



which on rewriting in terms of the Legendre effective action T A gives (as in ref. [I 



d D q dC~ 1 (q,A) 



(27T) 



D 



dA 



(q> -q) 



(2.2) 



Ta is defined by Ta^] + |<^.C _1 .<^ = — Wa[J] + J.y, where now = 6W A /5J is the 
classical field. In the limit £ — > 0, the dC~ x (q, A) / dA term restricts the momentum integral 
to the spherical shell q = A. Therefore any other terms in ( |2.2j ) containing C~ 1 (q, A) [or 
C _1 (A, A)], become ambiguous in this limit, since they in turn contain 6(0), which is ill- 
defined To properly resolve the sharp cutoff limit we must therefore isolate all such terms 
and treat these more carefully Jl^]. Thus, following ref. |U||, we separate from the two-point 
function the field independent full inverse propagator ^(p, A): 



(p, pO = 7 (p, A) (27r) n S(p + p') + I>] (p, p'; A) 



(2.3) 



so that f [0] = 0, and drop from both sides of (|2.2| ) the field independent vacuum energy 
contribution. The subtracted form of (12.21) can then be written 



1 

— tr 
2 



^(c- 1 +7 )- 2 .f. (i + [G-i+ 7 ]-i.r _1 



(2.4) 



We will assume that the square-bracketed [C _1 + 7] term, buffered as it is on both sides 
by r [<£>], 'almost never' carries a momentum p = A, i.e. that such points form a set of zero 
measure for the integrations in ( |2.4|) . (The validity of this 'zero measure assumption' will 
be addressed in sect. 5.) Then, from the definition of C _1 (p, A) we may safely take the 
limit 



lim 



1 



G(p,A) 



(2.5) 



£ -o C~ l (p, A) + 7 (p, A) 
where G(p, A) = 6(p — A)/ r y(p, A) is the I.R. cutoff full Green function. To complete the 
sharp cutoff limit of (|2.4| ) we need to note 

1 



,. dC l (q,A) r „_i, AN , ANl -2 d 



lim 



dA {e^o C - 1 (q,A)+' 1 (q,A) 
5(q-A)Mq,A) , 



9(0) is by no means simply given by 6(0) = \ [13] ! 
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where the last line results from taking the limit in a similar way to (|2.5|) . Substituting this 
and (|2.5|) in (|2.4| ) we obtain the sharp cutoff flow equation: 

d _ f f d D q 5{q - A) 



/(&S 1 [ f - (1 + G - frl ] ( "'- q) ' ™ 



<9A 2„ 

Under the field rescaling (|1.2[) we have that T i— > A 2 T, 7 i— > A 2 7 and G 1— » G/A 2 , as trivially 



follows from ( J2.3[ ). We see that the sharp cutoff flow equation ( |2.6| ) indeed has manifest 
invariance under field rescaling ( |1.2|) , as promised. This equation may now be expanded as 
a power series in the field ip to yield an infinite set of coupled non-linear flow equations for 
the (n-point) one particle irreducible Green functions T(pi, • • • , p n ; A). As usual, in mo- 
mentum space we factor out the momentum conserving 5-function [r(p, — p; A) = 7(p, A)]: 

(27r)^( Pl + --- + p n )r( Pl ,---, Pn ;A)= /7 A W S " ( 2 - 7 ) 

<V(pi) • • • 5ip(p n ) 

Our interest in using these equations is to investigate continuum limits, which are 
found by approaching a fixed point as A — > 0. This justifies dropping all reference to 
the overall cutoff Aq, which we do implicitly (explicitly in ref. flT3|), since all momenta 
can be restricted to be very much less than Aq, as is clear from the trivial bound- 
edness of the momentum integral in ( |2.6| ). At the fixed point there is only one mass 
scale, namely A, and therefore it is helpful to write all dimensionful quantities in terms 
of this. We need only observe that the scaling dimension of the field (p, at the 
fixed point, is generally anomalous: d v = ^(D — 2 + 77), where r\ is the anomalous 
scaling dimension. We could scale A out of the infinite set of flow equations for the 
r(pi, • • • , Pn! A), deducing their scaling dimensions from Q2.7j) , but it is more elegant to 
scale A directly out of the functional eqn. (|2.6|) . Therefore, similarly to refs.[|14"|,|8||, we write 
q h- Aq, y?(Aq) ^ A^-^qVC I/jC^A^x)] 1- C^M, l(^P, A) ^ A 2 ~^(p,t), 
f[C -1 A d ^](Ap,Ap';A) i-> A 2 - D -^f[( / 9](p,p / ;t) and t = ln(A /A). Here A is the 
energy scale at which the 'bare' Tq[ip] will be defined (Ao « Aq), and the number 
( = (47r) D / 4 v /r( J D/2) is chosen for convenience. (The physical meaning of To[tp] is dis- 
cussed in sect. 4). In addition we perform the radial part of the q integral in ( |2.6| ). The 
result is: 

(|+^ + A p )r4^] = -^([f.(l + G.f)- 1 ](q,-q)) . (2.8) 
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In here, = <p-j- is the 'phi-ness' counting operator: it counts the number of occurences 
of the field (p in a given vertex. A p may be expressed as 

r d D v d S 

^lp?"*V* (2 - 9) 

and is the momentum scale counting operator. The angular brackets stand for an average 
(i.e. normalised angular integration) over all directions of the unit vector q = 1. The full 
propagator G(p, t) is now given by 

G(p,t) = d(p-l)Mp,t) . (2.10) 

Using (1 + G.r) -1 = 1 — G.T + (G.f) 2 — ■ ■ • in ( |2.8|) and expanding in ip we obtain the 
infinite set of scaled flow equations for the n-point Green functions: 

d n d \ 1 

d- t + J2 p ^d p ^ + nd ^~ D J r(pi,---,Pn;0 = ^-^(r(q,-q,pi,...,p n ;t)> 
--fTT) 2 <r(q,-q-Pi,/i;t)G(|q + PiM)r(q-P 2 ,-q,/ 2 ;t)) 
+ ^ (r(q,-q-P 1 ,/ 1 ;t)G(|q + P 1 |,*)x 

T(q + Px, -q + P 2 , h- 1) G(|q - P 2 |, t) T(q - P 2 , -q, / 2 ; *)) + • • • . 

(2.11) 

These are illustrated in fig. 2. P^ = J2 Pk eii P fc ' an< ^ ^2{h h} h ■■■ i m * s a sum over a ^ 
disjoint subsets Ii PI Ij = such that Ul^i 1% = {Pi' ' " ■> Pn}- The symmetrization 

{ii, 7 2 } means this pair is counted only once i.e. {ii, I 2 } = {/2, Ii}- Evidently the expan- 
sion stops at the term where all vertices have their minimum number of legs, i.e. at the 
n th term in general, or the (n/2) th term in the (p «-> —<p invariant theory. Their structure 
is identical to that given in ref.[^3f, except for scaling out of A and the redefinitions to 
exhibit field reparametrization invariance (not discussed in ||13|1 ). The latter now appears, 
by (p.7|) ( p.l0| ), simply as 'n-point' invariance: 

r(pi, •••,?„;*) ^A n r( Pl , ...,p n ; t) n>2 . (2.12) 



Of course ( |2.8|) and (|2.11|) are still exact, however these are now in a particularly helpful 



form for considering the momentum expansion. 
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Fig. 2. The sharp cutoff flow equations for one particle irreducible vertices. Internal lines 
are full propagators. The black dot now represents restriction to momentum q = 1; the 
other propagators have an I.R. momentum cutoff so that p > 1. The symbol A represents 



the operator on the LHS (Left Hand Side) of eqn. (|2.11|) . 



3. The momentum scale expansion. 

The momentum expansion for the Green functions T(pi, ■ ■ • , p n ; t) is in terms of ho- 



mogeneous functions of non-negative integer degree[13 



r(pi,...,p n ;t)= J2 rim) (Pir--,Pn,t) 



m=0 



(3-1) 



such that r(™) (p Pl , • • • , pp n ; t) = p m T^ ( Pl , • • • , p n ; t) . 
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Or, before expanding in <p, we may write equivalently: 

oo 

r*M = £ r * (m) M where ^ m) M = ("» - ^)r! m) M (3.2) 



m=0 



[hence the name "momentum scale counting operator" ( |2.9|) . The extra factor of D arises 
from the momentum conserving 5-function in (|2.7|) .] All other external momentum de- 



pendence in (|2.11|) may be expanded as a series in integer degree homogeneous functions 



by introducing the "momentum scale" expansion parameter p, and expanding as a power 
series in p, e.g. the one particle irreducible Green functions generally are reexpanded in p 
via r( m )(q + pP, — q — pP', ■ • •), after which p may be reset to one. The I.R. cutoffs thus 
have expansions as: 

6>(|P + q| - 1) = 0(q.P + P/2) = 6>(q.P) + V -r (P/2) m 5^ m ~ (q.P) , (3.3) 

m=l 

where we have introduced the unit vector P = P/P (and used 6>(x — 1) = 9(x 2 — 1). 
<5^ m ^(x) means the m th derivative of the 5-function with respect to x.) Of course the 
fact that all external momentum dependence in (|2.11|) can be expanded in this way in 



terms of homogeneous functions of non-negative integer degree, guarantees that such an 
integral momentum scale expansion is self-consistent (as opposed to needing for example 
homogeneous functions of fractional degree). A systematic sequence of approximations - 
the 0(p M ) approximations - now results from replacing the sum Ylm=o * n dHUXlLl) by 
Sm=o> a sum U P t° some maximum power of momentum scale M = 0, 1, 2, • • • , substituting 
this into the flow equations (|2.£|) ( PTTT| ), and expanding the RHS (Right Hand Side) of these 
equations up to the same order in momentum scale. 

If the r(pi, • • • , p n ; t) were analytic functions of momenta at the origin = 0, then 
they would have a Taylor expansion in the momenta, equivalent to a derivative expansion 
in position space, and therefore odd m T^ m "> would vanish and even m r< m ) would be 
sums of Lorentz invariant products of m momenta. Such an expansion is not however 
possible for sharp cutoff P2]JT^| , a point we reemphasise here. The sharp cutoff induces 
non-analyticity at the origin of momentum space. This is already evident in the expansion 
( P-3| ): the odd m terms do not vanish and the 0(p )'th term, 0(q.P), is still a function of 
the angle between between P and q and not simply a constant as would be required by 
analyticity. Remarkably, this non-analytic dependence cancels out at 0(p°), as we will see 
in sect. 5. For an expansion in small momenta beyond O(p ), there is no such serendipity, 



10 



and the best that one can arrange is an expansion in momentum scale of integer degree 
as in ( |3. 1| ) ( pT2|) . This can be readily confirmed by computing the momentum expansion 
beyond lowest order, on perturbatively evaluated Green functions (with I.R. cutoff internal 
propagators) e.g. the 4-point function to one loop in theory (as was done in effect in 
ref. fl3[] and reproduced here in appendix A, for completeness). At this point we should 
also recall that this non-analyticity is only a technical problem, and does not indicate that 
there is anything fundamentally wrong with a sharp cutoff [|13| . 



4. Uniqueness and locality. 

So far we have indicated how to form the 0(p M ) approximations to the flow equations 
( |2.11|) . From their structure as first order differential equations in t, it is clear that they 
serve to determine uniquely the 0(p M ) approximation to the Green functions T^ m \ m = 
0, • ■ ■ , M, at 'time' t + St, from their values at time t, and hence uniquely from the 'initial' 
values of the r^ m ^ at time t = (i.e. A = Ao). It is here (without loss of generality) that 
we must recover the notion of locality. The point is that we want the full theory ( |2.1| ), 
in the limit A — > 0, not to suffer from these spurious non-analyticities at zero momenta 
(which correspond to spurious non-local behaviour in position space). In the continuum 
limit, this behaviour should be purely a consequence of inserting a sharp low energy I.R. 
cutoff in the path integral, and cancel out once the momentum modes with p < A are also 
included. Locality would usually be implemented by insisting that the bare action S\/ 
be local (that is, be expressible in terms of a derivative expansion). However this would 
require reintroducing the overall cutoff Ag, complicating the flow equations somewhat ||i~3H , 
and this we wish to avoid. Therefore we will impose instead, at A = Ao, that the 'bare' 
Legendre action To[(p] be local|T^]. By universality (of the continuum limits) we do not 
expect that this change matters for renormalized quantities, since it amounts to some 
minor alterations in how the U.V. cutoff is imposed. [Feynman diagrams vanish if any 
internal momentum is larger than Aq, while otherwise these diagrams vanish only if all 
internal momenta are larger than Ao ||13|| . The reason for this difference is simple: the 
latter is the complement to the requirement that is imposed by the infrared cutoff when 
A = Aq namely that the contribution would be non-vanishing only if all internal momenta 
are greater than Ao- That the weaker U.V. constraint at Ao already fully regularises the 
theory, is clearest from the boundedness of the momentum integrals in (|2.6| )( pT8| ) ( |2.11| ).l 
Using (|2.11| ) to evolve away from t = 0, non-analyticity will be generated by the terms in 
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however this non-analyticity is now fixed, up to the usual arbitrariness in the choice 
of local bare action, and acceptable in the sense that it can be expected to disappear as 
A//z — > (where \i <K Aq is some typical scale of low energy physics). This is a full 'in 
principle' solution to fixing the momentum dependence of the T^ m \ but, providing the 
angular averages on the RHS of ( |2.11| ) can be performed, we find in practice that much 
more can be done, as follows. 

One method is simply to ansatz the general momentum dependence of the (with 
appropriate t dependent coefficients). If the ansatz solves the equations ( |2.11|) , in the sense 
that the momentum dependence of both sides of the equation can be made to agree (within 
the 0(p M ) approximation) by appropriate choice of the flow equations for the coefficients, 
and if the ansatz is sufficiently general that it matches the local bare action [up to 0(p M )} 
at t = (by appropriate choice of boundary conditions for the coefficients^) , then it is 
the solution, by the uniqueness properties of the first order in t differential equation (as 
discussed above). 

Another method is to determine the solution systematically as follows. We have that 
at t = the are no n- vanishing only for m even, and the latter are generally arbitrary 
sums of products of m momenta [subject only to the general constraints of Lorentz in- 
variance, permutation invariance, momentum conservation etc. as follows from e.g. (|2.7|) 1. 
Substituting these general expressions for T = ^2 m= Q^ m ^ (P-Hl); an< ^ performing 



the momentum scale expansion, we can determine the general form of the non-analytic 
parts in r ( - m ) at some small time later t = St. From this we can deduce at the 'linearised 
level' the general 0(p M ) form of V at any time t by replacing those coefficients of order 
St by general coefficients. Now we substitute this expression into the RHS of (|2.11|) and 



hence deduce the general form of the 0(p M ) T at the 'quadratic level'. Iterating, we find 
in practice that this procedure converges and in this way we deduce the general 0(p M ) 
form of the r( m ) as a linear superposition of a finite set of 'basis' functions of momenta. 
Only the coefficients carry the t dependence and of course those multiplying non-analytic 
momentum functions, must vanish at t = 0. In the next two sections we give examples of 
0(p M ) approximations. 



10 Of course any terms which are then zero for all t, are dropped. 
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5. The local potential approximation. 

We will now work out the 0(p°) approximation. Since truncations of the field depen- 
dence tend to have limited accuracy and reliability |§, we make no further approximation. 

It is here that we must discuss the validity of the zero measure assumption, made 
below ( |2.4j) , which is equivalent to assuming that the ill-defined p = 1 propagator G(l,t) = 
9(0) t) [from (|2.10| )1 appears only at a set of points of zero measure in (p78|). [If this is 
the case then the result is the same whatever finite value one assigns to these #(0)'s, after 
performing the integrals in ( |2.8[ ).] Actually this assumption is plainly wrong when <£>(x) 
obtains a non-zero spatially independent vacuum expectation value v?( x ) = (v) HHL s i nce 
then T[ip] has a non-zero 'diagonal' part i.e. a part proportional to <5(p + p'). In ref. [13[ we 



handled this by redefining 7 in ( |2.3[ ) to absorb the diagonal part. One of the main results 
of this paper however will be to show that this procedure is not necessary. Instead, we will 
show that at 0(p ), the vanishing external momentum limit of the sharp cutoff equations 
(i.e. considering y?(x) — > (cp), after having taken e — > 0), is unique and equal to the result of 
taking the e — > limit of the equations where all momentum dependence is discarded first 



(as was done in refs.|T^,|S|,|T2||). In other words we will show that these two limits commute. 
As we will see, this fact is by no means trivial however: the agreement is obtained only at 
the end of the computation and in very different ways in the two methods. What is more, 
the 0(p ) approximation exactly coincides with the local potential approximation |4|] ( |1 . 1|) , 
providing further justification for the ubiquitous use of the latter PHl2| by demonstrating 
that it is simply the start of the systematic sequence of 0(p M ) approximations. Where 
derivations of the local potential approximation appear in refs. [fIHl2| > they are all of the 
same type as that in refs.|]I|^,^[: utilising a finite width momentum shell (with sharp 
edges) and discarding the momentum dependence first, before allowing the width to tend 
to zero. This is equivalent in effect to the derivation in refs.]T3|J^] if one specializes to 
a cutoff function that simply linearly interpolates between constant values outside the 
shell i.e. to 9 e (q, A) = (q — A)/e for \q — A| < e. The important point is that in all 
these derivations it is not possible even in principle to improve on the approximation of 
discarding all the momentum dependence, while the latter is done before letting the width 
e — > 0. This is simply because any non-zero external momentum p will satisfy p >>> e as 
e — » 0, but proper consideration of this regime precisely corresponds to interchanging the 
limits. Therefore any attempt to go beyond the simplest approximation of discarding all 
momentum dependence, must first address the question of whether this exchange of limits 
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still yields the same answer for the simplest approximation. As we have already stated, 
this turns out to be non-trivial but true. 

At 0(p ) the local bare action can only include interactions in a bare potential Vo(ip) 
which however is taken to be completely general to begin with.0 By either of the methods 
discussed in sect. 4, this suggests that we start by ansatzing 

VM = Jd D x{\{d^) 2 + V t ^)} . (5.1) 

This will turn out to be the complete solution to the momentum dependence at O(p ), 
i.e. the vertices T(pi, • • • , p n ; t) with n > 2, which are in principle still functions of the 
angles between the external momenta, may be taken to be simply constants. Note that 
the kinetic term will not get corrected at O(p ) [or 0(p 1 )} so carries no t dependence. The 
normalization for the kinetic term is conventional, but would require (and could at this 
order be given) separate justification here if field reparametrization invariance had been 
broken, since then the normalization in general affects the results [pO| , pT6| , jT^pTTjl . However 
in this case we have immediately that all normalizations are equivalent, by the invariance 
(|j~2|). Equation ( |5.1| ) implies that 

7 (p,t)=p 2 + K'(0) (5.2) 
and f[^](p,-p-P;t) = f(P) = y'rf D x[K'(^(x))-K'(0)]e' RX . (5.3) 

Here, prime refers to partial differentiation with respect to (p. The latter formula follows 
directly from 6 2 r/6<p(p)6<p(p') = jd D xd D y e ~^ P ' X + P ' y ^ 2 r/^(x)^(y). We have in- 
troduced the notation T(P) to emphasise that T depends only on the sum of its two 
momentum arguments; this is a consequence of dropping all momentum dependence in 
the vertices. Substituting this, ( |5TT| ) and (|2.10| ) into (pT8|), expanding the RHS of ( p78| ) via 

(1 + G.f ) _1 = 1 - G.f + (G.f) 2 , and dropping all but the O(p ) terms in 7 and 9 in 

the RHS (c.f. (|3.2| ), (|3.3| ) and the discussion inbetween) one obtains: 



\ Jd D x(d^ = 



11 Higher terms in momentum scale can be added, but at O(p ) there is no feedback into these 
terms from the RHS of the flow equations ( |2.8| ) ( |2.11| ) . The LHS of these equations then constrain 



these terms to be singular in <p (which is not allowed) or vanish. See ref. [14] or ref. ||15|| for further 
explanation. 
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implying r\ = for this approximation, as expected, and: 
J d D x { J^M + \{D- 2)<pVl{<p) - DV t (<p)} = 

f(o) ^ (-iy r d D P!d D p 2 ■ ■ ■ d°p r 



(5.4) 



Xt) ^<y(h 



7 (1, t) ^ 2 j(l,tyj (2n)^-D 



r(Pi)r(p 2 )---r(p T 



x (0(q.Pi) 0(q.[Pi + P 2 ]) • • • 0(q.[Pi + P 2 + • ■ • + P r -l]) > , 

where P r = — Y^k=i^*k- The r th contribution can be represented diagrammatically as 
in fig. 3, and alternatively follows from the O(p ) parts of contributions in fig. 2, after 
integration over the fields <p(pi) and resummation over n. 




q+ R +P 2 + - + P M 



Fig. 3. The r th contribution to eq.( |5.4j) , represented diagrammatically. Recall that the 
small dot indicates averaging over all directions of the propagators momentum q, which is 
resticted to be of unit norm. 

The above average over the unit vector q is indeed a non-trivial function of the angles 
between the momenta P/., but by utilising the cyclic symmetry of fig. 3 or (|5.4|) with 
respect to the external momenta, we can effectively replace its contribution by 1/r, as 
follows. We rewrite the RHS of ( |5.4| ) as: 



£(-ir 



r=l 



n ( t^7(M) ) {2lv)D6 [ti Pr I ( ' (q ' Q2) " (q - Q3) ' ' ' " (q - Qr: 



(5.5) 
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where we have introduced the partial sums Q fc = Ej=i p j for 2 < /c < r, Qi = 0, and 
defined the angular average to be unity for r = 1. Taking the average of the r integrals 
obtained by cyclically relabelling the integration variables Pi, • • • , P r (geometrically cor- 
responding to the r equivalent diagrams obtained from fig. 3 by placing the small dot on 
each of the r propagators) we obtain: 

u^^)^Mt p Mifi^^-QA)) . (5.6) 




* (2tO d 7(M) 



v fc = l / \j = l fc=l 



Now, generically for any given vector q, one of the q.Qfc, say q.Qi, will be less than all 
the others i.e. q.Qi < q.Qfc V7c 7^ i G (1, • • ■ , r). It follows that for this q, the j = i 
contribution to the above j-sum is unity, and the j 7^ i contributions vanish. Therefore, 
except for a set of points of zero measure for the q = 1 angular average (defined by 
those vectors q such that q.Qi = q.Qfc for some k), Y^j=i Tlk^j ^(<l-[Qfc — QjD = 1- 
Substituting this, and (|5.3|) , and writing (27r) D 5(^^ =1 P r ) = f d D x exp — i(x.Pi + • • • + 
x.P r ), we have immediately that equals - / d D x J^^Li 7 W"{ Q ) ~ v t"(v)Y h(l, tf ■ 
Using (|5.2[ ), this may be resummed to J d D x {in [l + V t "(ip)\ — ln7(l, t) }. Here, the second 
term, which corresponds to a vacuum energy, may be renormalised away by redefining the 
field independent part of the potential as V t ^ V t + £(t), where, by (|5.4j), we require 
88/ dt — D8 = — hi7(l,t). It is neater to notice however, that a 7 dependent term was 
discarded in ( |2.4| ) as a result of dropping from (|2.2|) the field independent vacuum energy 
contribution. If this 7 dependent term is kept, it is straightforward to show that it precisely 
cancels the — J d D x In 7(1, t) above. Hence finally, we have from (|5.4|): 



+ ~ (D - 2) <pV/(<p)-DV t (<p) = In [I + V/' &)] . (5.7) 

This coincides with the local potential approximation0-|T^] to the Wegner-Houghton 
equations [fflj. 



6. Beyond lowest order. 

We now proceed to the 0{p 1 ) approximation. Following the iterative method described 
in sect. 4 we try to proceed by initially ansatzing a momentum independent form for Ft as 
in ( |5.1| ), (|5.2| ) and (|5.3| ). We then substitute these expressions into (|2.8| ) to see what form 
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the 0(p 1 ) T t takes at the 'linearised level'. We see from expanding j(p,t) in (|2.10| ) and 
using (|3.3|), that amongst other terms a set of averages of q.Qfc (k = 2, - • • ,r) of the form 



(0(q.Q 2 ) • • •0(q.Q fc )q.Q fc 0(q.Q fc+ i) • • -0(q.Q r ) > , (6.1) 

have to be computed [where we have used the notation introduced in (|5.5|) ]. While we 
were able to replace the analagous expression in (|5.5| ) by a constant, using (cyclic) sym- 
metry arguments, we have found the permutation symmetry of the diagrams fig. 3 to be 
insufficient, when r > 3, to allow a simple replacement for the above expression. Of course 
the above expression can be evaluated exactly, in principle, but the result will be a very 
complicated 0(p 1 ) function, so that further insight is necessary to make progress at this 
level. 

For this reason we now restrict the discussion of momentum dependence to cases where 
r < 3. This may be 'justified' by approximating the flow equations by truncations of 
the field dependence |13[ of the O(p m>0 ) terms to six-point functions or less (in a ip «-> 
—f invariant theory), i.e. discarding ~ p 8 terms and higher in non-zero momentum 
scale pieces. However, we already know that truncations of the field dependence yield 
results that do not converge (beyond a certain order) and generally also result in spurious 
fixed points M. Nevertheless, we will discuss these truncations in order to provide explicit 
examples of the general considerations of the previous sections. As we will see, the results 
are not impressive. In particular they are much worse than the results of the momentum 
expansion with smooth cutoff in ref.[^4f. We attribute this to the truncations of the field 



dependence, because this was not employed in ref.|14[], the sharp cutoff momentum scale 
expansion results in a rapidly convergent (numerical) series in a model two-loop calculation 
( fll3|1 and app.B), and intuitively one would expect a better convergence for sharp cutoff 
than for smooth [|T3|] . Although we cannot at this stage rule out the possibility that the 
sharp cutoff momentum scale expansion is 'at fault' at the non-perturbative level, we 
suspect that truncations of the field dependence of the kinetic term K(tp) in the smooth 
inj would produce similarly unimpressive results, which would then strongly support 
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our contention that it is also only the truncations that are at fault here. 

Because we will be working to 0(p m ) with m > 0, only for certain n-point functions, 
we must generalise the discussion of previous sections to the cases where m can depend 
on n, i.e. where the equations ( |2.11|) are expanded to different order in momentum scale 
depending on the value of n. This generalisation is however straightforward and will be 
explained within the examples below. 
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We concentrate again [ TJ ] || on the critical exponents for the non-perturbative Wilson 
fixed point in three dimensions: the results will provide a simple measure of how good 
these approximations are, in a non-perturbative setting. In three dimensions, averages of 
functions / that depend on only one component of q, say q.P as in (|3.3|), satisfy: 



</(q.P)> = i fdzf(z) . (6.2) 

It is of course straightforward to give the analogous formula in any other dimension, four 
for example [[TBI, and all averages we consider can be reduced to this (or its analogue) by 
symmetry arguments, in any dimension. 

Consider first the simplest possible truncation that will allow non-trivial momentum 
dependence: discarding all but the two-point and four-point Green functionsil! By ( |2 . 1 1| ) 
this gives 

(at + p d~ P ~ 2 + t] = WT) <r(q ' -q ' p ' ~ p; ^ (6 ' 3a) 

d 4 d 

+ XX ^ + 1 ~ 2? ?) r (Pi> P2 5 P3, P4; *) = (6.36) 



i=l 

2 



7(1.*) 



r(q, -r 12 , Pl , p 2 ; t) G(r 12 , t) r(r 12 , - q, p 3 , p 4 ; t) + (2 <-> 3) + (2 <-> 4)) 



where = q + Pi + Pj , and in the second average the last two terms are the same as 
the first term but with indices swopped as indicated. We have thrown away the six point 
contribution ~ (T(q, — q, pi, p 2 , P3, p 4 ; t)) to ( |6.3fc| ). Working to 0(p 1 ) we start again by 
ansatzing first that the four-point function is momentum independent. The analogue of 
( |Q| ) is now, using (|6.2|), simply ^(q.P)q.P) = P/4, where P = pi + p 2 , pi + P3 or 
pi + p4 depending on the term on the RHS of (|6.36|) . Similarly the 0(p 1 ) term in ( |3~^ ) 
produces these P's. Therefore at 'linearised level' we have that 

r(pi, P 2,P3, P4;*) = O!o(*) + Oil(*){|Pl + P2| + |Pl +P3| + |Pl +P4|} , (6-4) 

where flow equations for the are to be determined. Substituting this into eqn.( |6.3a| ) 
results in p dependence of the form (|q + p|) = l-|-]9 2 /3 [by fl6.2|) ], therefore a quadratic 



12 This is closely similar to a model computation in ref. |]T^| . Were we also to throw away all 
momentum dependence [H it would be a differential version of the "5 4 model" [E2[. 
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ansatz remains sufficient for ^(p, t) at this order. Substituting (|6.4| ) back into ( |6.36|) results 
in contributions e.g. of the form 



r(q + Pi +P2,-q,P3,P4;£) = «o +ai{2+ |pi + p 2 | +q.(pi + P2)} + 0(p 2 ) , (6.5) 

where we have used momentum conservation and q = 1, (and similarly for the other such 
contributions). Therefore the RHS of (|6.35|) reproduces the general form of (|6.4|) ; our 
iterative procedure has closed at the linearised level. To obtain an approximation for r\ we 
will work to 0(p 2 ) with j(p, t), writing 

7 (p,t) = a(t)p 2 + a(t) , (6.6) 

and expanding (|6.3a|) to 0(p 2 ), but keep to 0(p 1 ) with T(pi, • • • , P4; t). Substituting 
and (|)]6|) into ( |6.3|) we obtain: 



dt + {v ~ 2)a = "7+7" (6 - 7g) 

dt 3 a + a 

da /a + 2ai\ 2 

—- + (2r7- 1 "0 = -3 6.7c 

dt \ a + a J 

dcti n 1 (a + 2a 1 )(aa - aa + 8cxai + AaaA (n „ A 

— -+2r/a 1 = -- . 6.7d 

at 2 (a + a) 6 

By the locality requirements of sect. 4 we have that the bare cr(0), a(0) and «o(0) are a priori 
arbitrary, but «i(0) = 0. By the field reparametrization invariance ( |1.2|) , equivalently 
n-point invariance (p,12| ), we can impose the conventional wavefunction renormalization 
condition a(oo) = 1. (The full renormalized quantities are given by their t = 00 values.) 

Had field reparametrization invariance been broken (as happens generally under such 
approximations [pO| , |il||Tl~f ) , one would discover the following problem: the approximate 



physical quantities would differ for different values of a(oo), even though they do not in 
the exact formulation [as follows by using (|1.2|) 1. This problem only potentially arises 
when we seek to go beyond the local potential approximation, and is not yet generally 
appreciated in the corresponding literature, largely because a(oo) (or its analogue) is set 
to unity anyway (in effect by the tuning of bare parameters or the form of the cutoff), 
even when the justification for this specialization is missing. To date, only two forms of 
cutoff function are known that allow a field reparametrization invariance after the mo- 
mentum expansion approximation: the sharp cutoff studied here, and power-law smooth 
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cutoffs 0|15|,|TTp C(q,A) ~ (q/A) 2K M The problem is particularly clear if we study the 
Wilson fixed point p2[. We assume that choices can be found for a(0), a(0) and ao(0) that 
lie on the critical surface. Then as t — > oo, the four quantities cr(t), a(t), ceo(t) and cti(t) be- 
come t independent and eqns.(|6T7|) collapse to four simultaneous algebraic equations. The 
problem is however, that with 77, there are five quantities to be determined. Without field 
reparametrization invariance, this allows r\ to be determined only as a function of 0(00) 
(say). With field reparametrization invariance, these equations correspond to non-linear 
eigenvalue equations for the field's anomalous dimension 77, which therefore (generically) 
can have only certain discrete values - as it should be. 

Solving the eqns.(|677|) at a fixed point a{t) = a, a(t) = a, etc, where they reduce to 
algebraic equations, we obtain that ao, a.\ and 77 can be written as polynomials in a, which 
itself satisfies the quintic 

3a 5 + 92a 4 - 154a 3 - 159a 2 + 74a + 8 = , 

except for the Gaussian solution a = oti = r\ = 0. By the scaling (|1.2|)( |2.12"D , we have 
set (in all cases) a = 1. The quintic has five real solutions a = —32.2, —0.953, —0.0917, 
0.441 and 2.15, however only a = —0.0917 gives a positive cto (and hence a physically 
stable potential). We assume that the domains of attraction of the other fixed points do 
not include physically sensible bare actions with a(0) > 0, «o(0) > and cti(0) = 0.0 
Using a = —.0917 we obtain ao = .103, a.\ = .0307 and 77 = .0225. This last value, 
being universal, can be compared to the best determinations: r\ = .032(3) or .038(3), 
where the first derives from perturbation theory in fixed dimension and the second from 
e expansions |2^]. Computing the linearised perturbations of ( |6.7|) about this fixed point, 
we find a numerical matrix eigenvalue equation^! whose eigenvalues are A = 1.88, 0, 
and —.390 ± .630i The most positive eigenvalue yields v = 1/A = .532 by standard 
arguments]^], differing little from the S 4 model[]22| {y = .527), and should be compared 
to the best determinations v = .631(2) |[25|| . The zero eigenvalue corresponds to the 



13 It can be shown that for other choices the field reparametrizations [B4[ only leave the flow 
equations invariant if they can carry non-trivial momentum dependence, and this is then broken 
by momentum expansion to finite order. 

14 In general however, there are no reliable ways to reject the spurious answers resulting from 
truncation of the field dependence [||. We note that in this case, the other values of a also result 
in wildly large values for the other parameters. 

15 N.B. a(t) = 1 + 5a(t) where 8a(t) is not fixed. 
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reparametrization symmetry ( |1.2|) and is therefore redundant [ 23 , H] . Finally, the irrelevant 
complex eigenvalues represent the (first) corrections to scaling for this truncation ui = 
—X = .390 ± .630z, which should be compared to the best determinations u = .80(4) pSfl . 
This eigenvalue is the most affected by the approximation, as should be suspected. These 
results could probably be improved by making the truncation around the minimum of the 
fixed point effective potential [|17|,|9[|. Our point here however, is to begin with the simplest 
possible non-trivial example of the general method of momentum scale expansions. 

We turn next to more sophisticated truncations, although we are limited to truncations 
of the field dependence in the momentum dependent sector, as we pointed out above. First 
let us generalise the above example by avoiding the truncation of the potential, that is we 
include all the higher n-point functions T(pi, ■ ■ ■ ,p n ;t) with n > 6, at O(p ), with the 
corresponding n > 6 flow equations ( |2.11|) also evaluated to O(p ). In this case it is clear 
that the ansatz that these higher point functions are independent of momenta, must satisfy 
these equations, since the equations differ from those of sect. 5 only in the following ways: 
a(t) 7^ 1, 7] 7^ 0, and the four-point function appears in O(p ) expressions as ao(t) + 2ai(t), 
as follows from ( |6.5|) (and its generalisations). By sect. 4, this is then the unique answer 
for this truncation, up to the choice of bare potential Vo((p), and bare kinetic term a(0), 
for the local bare action. The resulting flow equations are the higher momentum scale 
equations (|6.7fc, d|) and the flow for the potential which may be resummed to: 

d - Vtlv) + 1(1 +V) <pV{[<p) - 3V t (<p) = In [o(t) + «i(% 2 + V?(<p)] . (6. 



dt ^ 2 

The changes between ( |5.7| ) and this equation just follow from the differences with sect. 5, 
listed above. At the fixed point we again set a(t) = a = 1, using (|1.2|). Equations ( |6.7a, 5, clj ) 
are unmodified by the inclusion of the O(p ) n > 6-point functions [( |6.7a|) being subsumed 
in ( 6.8a )1, and may be solved in terms of a as «o = (1 + cr)([a — 3]rj — 2a), a.\ = |(1 + a)r] 
and 

9 + 7a + 2a 2 ± ^(1 + a)V35 - 8a .„ c ,. 

77 = a 5 • 6.86 

' 6 + 21a + 14a 2 + a 3 V ; 

At a fixed point, Vt((f) = V((p) therefore depends only on the choice of a which may be 

regarded as the <p = initial condition, together with the evenness constraint V(0) = 0. 

As explained previously Q one finds that all but a discrete set of solutions of ( |6.8| ) are 

singular (at some finite (p) and therefore unacceptable. We find only two non-singular 

solutions, being the trivial Gaussian solution V(<p) = 0, and an approximation to the 

Wilson fixed point. Quantitative methods for computing the latter are described in ref.H. 
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(See also refs.|n],|15]. We used the second of the better expansion methods in 0.) It is 
significant that all spurious results have disappeared: this is a benefit of not truncating the 
field dependence in the potential We find a = —.214 and r\ = .0660. Solving for small 
perturbations about this fixed point we deduce v = .612 and u = .91. As expected, 
this is a significant improvement over the crude truncation taken previously. These results 
are an improvement over those of the pure local potential approximation ||[7]||]: rj = 0, 
v = .690 and u = .595, but are not as good as the 0(p 2 ) smooth cutoff approximation [|l~4| . 
Again, we attribute this to the effects here of truncation of the field dependence in the 
0(p 1 ) part. 

Now we try to improve the results by truncating the 0(p 1 ) part at the maximum prac- 
tical value of n = 6. Here the computation starts to become quite involved, but once the 
general form of the 0{p 1 ) dependence is determined, it proceeds similarly to the previous 
example. Since our primary purpose is to illustrate the techniques of the momentum scale 
expansion, we indicate here only how to determine the momentum dependence, and then 
present the results. Since we now wish to include the 0(p 1 ) dependence of the six-point 
function, we need to work explicitly with the n = 6 case of the flow equations ( |2.11| ). 
We can determine the 0(p 1 ) dependence by iteration, substituting the previous solutions 
for the two-point ( |6.6j ) and four-point functions (|6.4|) . The six-point function's momen- 
tum dependence then arises in particular from the sum over all permutations of the legs 
Pi, • • • , P6 in the graph fig. 4, that is from the angular average: 

<r(q, -q - Pi, Pi, p 2 ; t) G(|q + PiU)x 
P(Pi>— >Pe) 

T(q + P l5 -q+ P 2 , p 5 , p 6 ; t) G(|q - P 2 |,*)r(q - P 2 , -q, p 3 , p 4 ; t)) . 

(6.9) 

Here P(pi, • • • , Pe) is a permutation of the six momenta Pi, • • • , P6- Pi = Pi + P2, P2 = 
p 3 + p 4 and P 3 = p 5 + p 6 . 

Expanding ( |6.9|) to 0{p x ) we find that, as well as averages that directly follow from 
( |6.2|) or sect. 5, the following averages need to be determined: 

P a (0(q.Pi)<y(-q.P 2 )) (6.10a) 
P 1 ^(-q.P 2 )o(q.P 1 )) (6.106) 
(0(q.P 1 )0(-q.P 2 )q.(P 1 -P 2 )) , (6.10c) 

where it is always to be understood that these contributions appear summed over the 
permutations "P(pi, • • • , Pe)- First note that ( |6.106| ) contributes the same as ( |6.10a| ), as 
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P 4 A 



Fig. 4. One graphical contribution to the flow of the six-point function. 

can be seen by substituting q i— > — q in the average and then permuting Pi and P2. 
(This permutation simply reorganizes the terms in the sum over permutations.) Now, 
substituting q 1— > — q in (|6.1Ua| ) gives the equivalent expression P 2 (6>(— q.Pi) S(— q.P 2 )), 
on using the evenness of the 5-function. Taking the average of this and (|6.10a|) , and using 
0(q.Pi) +0(-q.Pi) = 1, we deduce that (|6.10a|) equals ±P 2 (5(q.P 2 )) = This solves 

for ( |6.10a, 6f) . Using 6*(— q.P 2 ) = 1 — 6>(q.P 2 ) in (|6.10d ) we obtain the equivalent expression: 



(0(q.Pi) q.Pi) - <^(q.Pi) q.P 2 > - (<?(q.Pi) 0(q.P 2 ) q.(Pi 



The last term in this expression yields zero by the Pi «-> P 2 symmetry of the sum over 
permutations. In the second term we may replace P 2 by P3, using permutation symmetry, 
and hence replace the second term by — ^(6>(q.Pi) q.(P 2 + P3)). But this is just half the 
first term, by momentum conservation. Lastly, the first term is Pi/ 4, by ( |6.2|) . 

Computing the momentum dependence in this way we find that the six-point function 
is, to linearised level: 



r(pi,---,Pe;*) =g (t)+gi(t) |Pi + Pil 

pairs{i,j} 



(6.11) 
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where the flow equations for the gi are to be determined. Substituting this, (|6.4| ), (|6.6|), 
and constants for all 0(p ) higher n-point functions, we find the equations again close 
at the linearised level. The O(p ) parts resum into a partial differential equation for the 
potential Vt((p). The bare locality conditions require now «i(0) = gi (0) = 0. Determining 
the first few couplings at the Wilson fixed point, we find that rj can be expressed as the 
solution of a cubic with coefficients that depend on a, after which non-singular solutions 
of the fixed point equation for V (<£>), and perturbations around it, can be studied. We 
findS rj = .0591, v = .604 and 10 = .540 ± A7i. This estimate for ry has improved slightly 
compared with the previous truncation, but the estimates for v and u are actually worse. 

Finally, let us mention that we have extended the approximation to 0(p 2 ) for the 
four-point function. The only new average that presents itself is a term of the form 
(6 (q.P) qPq") = 6^/6, which may be evaluated in a similar way to that of ( |6.106| ). 
We find, by iteration as before, that the four-point function ( |6.4j ) aquires the new term 
d2(t) {pi +P3 + P4 } , otherwise all n-point functions have the same form as before. r\ 
is now given implicitly as a function of a through solutions of a complicated sextic. Unfor- 
tunately the sextic seems to have no sensible solutions (i.e. such that a < and rj > 0)1 If 
we drop instead to an O(p ) ansatz for the six-point function [i.e. set gi = in (|6.11D ], we 
find that the polynomial bounds the anomalous dimension to be r\ < 0.016, and therefore 
we must obtain worse results than our previous approximations. We again attribute these 
negative results to the, in general, poor behaviour of truncations of the field dependence 
- although, as discussed above, more research is necessary to confirm this conjecture. 



7. Conclusions. 

We have shown how the renormalization group for the effective action with a sharp 
momentum cutoff [0 , may be organised by expanding one-particle irreducible parts in terms 
of homogeneous functions of momenta of integer degree i.e. by performing a momentum 
scale expansion of the I.R. cutoff Legendre effective action. A systematic sequence, M = 
0, 1,2, of 0(p M ) approximations then follow. Significantly, the O(p ) approximation 
coincides with the local potential approximation^] and provides further justification for 
the ubiquitous use of the latter P-4T2| by demonstrating how it might be systematically 
improved. 

16 To simplify matters, we restricted our search only to the promising region of a where a < 
and rj > 0. 
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Note that the momentum scale expansion, although it is unambiguous to perform, is 
not an expansion in a small parameter: although the external momenta are regarded, for 
the purposes of the expansion, as small compared to the cutoff A, they are in fact set to 
values p ~ A when the results of the expansion are substituted into the RHS of the flow 
equations. It is a trivial but important point to note that this implies that the expansion is 
a numerical one which either converges or does not:0 in other words there is no hinterland 
of asymptotic convergence, where all of quantum field theory's small parameter expansions 
lie. It is not yet known in general whether the momentum expansion actually converges 
(although low order results compare favourably with other methods when such methods 
are available |19[|). We might expect that it generally does, because it corresponds to an 
expansion in an appropriately defined 'localness' of the effective Lagrangian, the expansion 
in momentum scale p corresponding in position space (for sharp cutoff) to an expansion in 
inverse powers of the average relative distance r between any two points in a vertex. Indeed, 
if there are large non-local contributions to the effective Lagrangian, then the description in 
terms of the chosen fields is probably itself inappropriate and indicates that other degrees 
of freedom should be introduced. This is very different from the cases where truncations 
of the field dependence itself could be considered appropriate: such truncations will only 
be valid if the field amplitude fluctuations can be considered small - but this is only true 
when mean field theory is a good approximation, i.e. precisely the regime where (weak 
coupling) perturbation theory is valid. In this way, we can understand why truncations of 
the field dependence in truly non-perturbative situations are somewhat unreliable^, while 
momentum expansions of the effective Lagrangian (with no other approximation) seem to 



be so successful [fl9f| . 

Another trivial but important point is to note that universal quantities, such as critical 
exponents, depend, at finite order in the momentum expansion, on the 'shape' of cutoff 
taken - whether it be sharp or smooth of some form C(p, A). More crucially, the sharp 



cutoff and power-law smooth cutoffs |14| are the only forms of cutoff known that preserve a 
field reparametrization invariance of the flow equations after approximation by momentum 
expansion. This is necessary for obtaining, with a given form of cutoff, a unique value@ 
for the fields anomalous dimension. 



17 Since the equations to be solved are non-linear, it is also possible that solutions fail to exist 
beyond a certain M. 
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i.e. universal, in the sense of being independent of the details of the bare action. 
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Note that, while Taylor expansions in momenta (a.k.a. derivative expansions) are 
possible for any finite width e of the (smooth) cutoff function, as e — > the dominant 
terms in these expansions are positive powers of ratios p/e, and diverge. The alternative 
expansion in momentum scale must be used in the sharp cutoff limit. As e — > 0, the flow 
equations become independent of the shape of cutoff and considerably simplify compared 
to their smooth cutoff cousins. As discussed in appendix B, the momentum scale expansion 
also appears to be rapidly convergent. There is a price to pay however: the homogeneous 
functions of momenta, which appear at each order of the momentum scale expansion, carry 
in principle much more information than Taylor expansions. Therefore it seems that one 
'pays' for the better convergence by absorbing more complication into the expansion at 
each order. Indeed, we did not succeed here in computing beyond O(p ) without truncating 
the field dependence. In the smooth cutoff case, it is clear that the expansion of the flow 
equations can be organised in such a way as to give differential equations only [ TA ] but here 
there is not such a simple action on the coefficient functions of the field. It is quite possible 
for sharp cutoff, beyond O(p ), that the resulting flow equations, if the field dependence 
is not truncated, would only be castable as integro-differential equations. One also needs 
(in principle) to check that a given fixed point behaviour can be reached from a local 
bare action: otherwise the behaviour may be the result of non-local physics that does not 
disappear as the I.R. cutoff becomes much less than typical low energy scales. 

Finally, we noted that the results of the 0(p°) approximation are well-defined for 
exceptional momenta, i.e. the limit — > is well-defined for any set of exceptional 
momenta - despite the fact that simply putting = results in 0(0) 's which are ill- 
defined. In our examples beyond 0(p ), in sect. 6, the limits are again well-defined. Let 
us mention that in all these cases, the results obtained are the same as would be obtained 
by taking the exceptional case first and then letting e — > 0. We do not have a proof that 
these properties are enjoyed to all orders in p n . If an ambigous result was obtained at some 
0(p n ), in the sense that the limit for exceptional momenta depended on the direction in 
momentum space in which this limit was taken, then one could presumably proceed by 
defining these quantities (if necessary) by the value obtained when the limit e — > is taken 
after taking the exceptional case. 



Acknowledgements 

It is a pleasure to thank the SERC /PPARC for financial support through an Advanced 
Fellowship. 



26 



Appendix A. Momentum scale expansion to one loop. 



In these two appendices we review some aspects of model results obtained in ref.|T3 
We rearrange them a little, to suit the present purposes. Working in four dimensions 
and with massless^l (<p <-> — ip invariant) A<^> 4 theory, one obtains to one-loop that the 
one-particle irreducible four-point function is given by 



i=2 



X 2V- / 3 d 9 fl(|q + Pi|-g) . . 



where Ao is the bare coupling, at tree level given by Ta [ip] = Jd'h; ^(d(p) 2 + Ao<^ 4 /4!. The 
momentum integral is restricted to A < q < Ao, as indicated. Pi = pi +Pi, and pi, • ■ ■ , P4 
are the four external momenta. The above expression corresponds to the standard s, t and u 
channel one-loop contributions. It may be obtained from one iteration of the flow eqn. ( |2.6|) , 
and takes the form expected - with a sharp I.R. cutoff multiplying each propagator - after 
a little rearrangement (c.f. sect. 4 and ref. |13|]). The momentum expansion of the flow 
equations corresponds, at this order in pertubation theory, to expanding (|A.1|) directly 
in small P^. We immediately see that a Taylor expansion is not possible, as a result of 
the function, whose expansion may be taken from ( |3.3|) . In fact using this, and Taylor 
expanding the denominator in ( A.l ), one obtains: 



The essential point is that this is an expansion in momentum scale Pi, and not a Taylor 
expansion in momentum components P?. In fact, the even powers of momentum scale are 
missing, in this example, to all higher orders. 



Appendix B. Convergence in a two loop example. 

In this appendix we go on to investigate the convergence of the momentum scale ex- 
pansion in perturbation theory. It is at two loops where one first has to confront the fact 
that the momentum scale expansion is not an expansion in a small parameter. In fact the 

19 This is the less favourable but crucial case: the massive case has better convergence properties. 
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Fig. 5. The two- loop four-point graph whose momentum expansion does not trivially 
converge. 



only contribution to the four-point function that does not trivially converge, is shown in 
fig. 5. 

It arises by iterating (|2.6| ) twice, and corresponds to the sum of two contributions, 
one where (say) the t channel contribution from ( |A.1| ) has two of its legs tied together by 
replacing Pi and A by p, weighting by the two outer propagators i.e. by ~ 1/p 4 , and then 
integrating over A < p < Aq, and another similar contribution constructed from the one- 
loop six-point function. Since the two contributions behave similarly under momentum 
expansion, we only investigate the former. Furthermore, we restrict our attention to the 
divergent parts, and can therefore set all external momenta to zero. (The finite parts in 
fact converge somewhat faster.) Therefore we are left to investigate a contribution of the 
form 

" A ° d 4 p 1 f A ° d 4 q 0(|q + p| -q) 



6A °A [2*y p*J p {2*y ^(q + P) 2 ' (R1) 

In the momentum scale expansion we are to expand the inner one-loop integral as we did in 
appendix A. Therefore we can simply substitute the result corresponding to ( |A.2| ). After 
performing the p integral one obtains the exact leading divergence as expected |fP3f, while 
the subleading divergence is given as a rapidly convergent numerical series: 

_/32 (4^ ln T' 

where fa = ^(8 + + + ' ' ') * s a contribution to the two-loop /3-function coef- 
ficient. The ratios of the partial sums to the exact answer fa = 2.568818, are r = 
.99130, .99956, .99996, • • •, corresponding to the 0(p M ) contributions where M = 1, 3, 5, • • •, 
and converge to 3sf already at (9(p 3 ), after which approximately an extra decimal place 
in accuracy is added with each new term. We expect that in this way, the sharp cutoff 
momentum scale expansion can be shown to converge at any desired order of the loop 
expansion. 
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Now we briefly compare with smooth cutoffs. The convergence of the momentum 



expansion with power law smooth cutoffs is rather subtle |fl4]| , in that the polynomial cor- 
rections to the full inverse propagator play a crucial role in the convergence of the mo- 
mentum integrals. Their momentum expansion convergence properties will be discussed 
elsewhere. Here we consider only a simpler example of an exponential infrared cutoff 
1/g 2 t— > 8(q, A)/g 2 , where 6(q, A) = 1 — exp(— g 2 /A 2 ), which can directly be compared with 
the above example. We might expect at least as good a convergence of the momentum 
expansion using this cutoff as with power law cutoffs. The smooth cutoff equivalent of 
O is 



where the one-loop integral I q is given by 



T _ fd 4 q g(|q + p|,A 2 ) d 

i ^A2^ 2 ^ 1 ( p 2 

1 — exp 



A 2 V V ) I V 2A 



The momentum scale expansion of I q results in convergent momentum integrals and gives 
of course the second line expanded as a power series in p 2 / A| . Performing the momentum 
expansion and the remaining integrals, one obtains the following closed expression for the 
corresponding @2- 



^n(n + l)\2j \2j 



^ n(n+l) V 2 

The partial sums up to n = M correspond to the 0(p 2M ) approximations. It is clear that 
this sequence converges. /?2 is different from its sharp cutoff value because, as described 
above, we have not included the term that follows from the one-loop six point function. 
We can compare ratios (of the partial sums to the full sum) however. In this case, these 
are r = 1.18269, .95273, 1.01432, • • •, corresponding to the 0(p M ) contributions where M = 
2, 4, 6, • • •. We see that convergence is about twice as slow as for the sharp cutoff momentum 
scale expansion, in this example. Slower convergence for smooth cutoffs is expected, on 



general grounds [13 
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